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Abstract 

We demonstrate that, for the case of quasi-equipartition between the velocity and the magnetic field, the 
Lagrangian-averaged magnetohydrodynamics a— model (LAMHD) reproduces well both the large-scale 
and small-scale properties of turbulent flows; in particular, it displays no increased (super-filter) bottleneck 
effect with its ensuing enhanced energy spectrum at the onset of the sub-filter-scales. This is in contrast 
to the case of the neutral fluid in which the Lagrangian-averaged Navier-Stokes a— model is somewhat 
limited in its applications because of the formation of spatial regions with no internal degrees of freedom 
and subsequent contamination of super-filter-scale spectral properties. We argue that, as the Lorentz force 
breaks the conservation of circulation and enables spectrally non-local energy transfer (associated to Alfven 
waves), it is responsible for the absence of a viscous bottleneck in MHD, as compared to the fluid case. 
As LAMHD preserves Alfven waves and the circulation properties of MHD, there is also no (super-filter) 
bottleneck found in LAMHD, making this method capable of large reductions in required numerical degrees 
of freedom; specifically, we find a reduction factor of 200 when compared to a direct numerical simulation 
on a large grid of 1536 3 points at the same Reynolds number. 

PACS numbers: 47.27.ep; 52.30.Cv; 95.30.Qd; 47.27.E- 
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I. INTRODUCTION 



When large-scale numerical simulations of astrophysical or geophysical magnetohydrodynam- 
ics (MHD) are desired, all dynamical scales of the physical system are rarely, if ever, resolved. For 
this reason, sub-grid- scale (SGS) modeling of MHD dynamics in the context of computations in 
the geophysical and astrophysical context is required. This modeling can be achieved implicitly, 
in the simplest example by employing a dissipative numerical scheme, or it can be done explic- 
itly by creating a Large Eddy Simulation (LES-see OJ] for a recent review). Explicit methods for 
MHD are not as pervasive as they are in engineering, or for geophysical and atmospheric flows. 
In fact, modeling for MHD is a relatively new field (see [Q, [2I). One problem with extending 
the LES methodology for hydrodynamic turbulence to MHD is that most LES are based upon 
eddy-viscosity concepts [1], which can be related to a known power law of the energy spectrum 



|4j] (although generalizations can be devised, see e.g. [50), or upon self-similarity. For MHD, the 
underlying assumption of locality of interactions in Fourier space is not necessarily valid [a LZQ 
(a contradiction of self-similarity) and spectral eddy-viscosity concepts [8] cannot be applied in 
a straightforward manner as neither kinetic nor magnetic energy is a conserved quantity and the 
general expression of the ener gy s pectrum is not known at this time u. 
Purely dissipative models |Jil 17] are inadequate as they ignore the exchange of energy at sub- 
filter scales between the velocity and magnetic fields and such models have been shown to suppress 
small-scale dynamo action IT 1 311 and any inverse cascade from the sub-filter scales II 1911 . A satis- 
factory LES for MHD has been proposed for the case starting with some degree of alignment 



between the velocity and magnetic fields 
to low magnetic Reynolds number 1121 , 







22, 



20|] • Other restricted-case MHD-LES are applicable 



23[|. Extensions of spectral models to MHD based 



on two-point closure formulations of the dynamical equations proposed recently look promising 
in the analysis of turbulent flows and of the dynamo mechanism [5]. Finally, though technically 
not an LES, there are also hyper-resistive models for MHD which require rescaling of the length 
(wavenumber) scales to a known direct numerical simulation (DNS) Q 1 8Q ■ 

One model which can be written as an LES is the Lagrangian-averaged MHD (LAMHD) equa- 



tions D24j, 



25L 



260 ■ It has been shown to reproduce a number of features of DNS . In two dimensions 



(2D) for Taylor Reynolds numbers (R\) up to ~ 5000 it has been shown to reproduce selective 
decay, the inverse cascade of mean- square vector potential, and dynamic alignment between the 
velocity and magnetic fields [27] as well as the statistics of small-scale cancellation [28] and inter- 
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mittency I29L In three dimensions (3D) at Reynolds numbers (Re) of ~ 500, LAMHD reproduced 
the inverse cascade of magnetic helicity (associated with the development of force-free magnetic 
field) and the helical dynamo effect 113011 - It has also been tested (up to kinetic Re ~ 3000, mag- 
netic Re ~ 300) for its ability to predict the critical magnetic Reynolds number for a non-helical 
dynamo at low magnetic Prandtl number Dili . LAMHD performed well in all these tests. Its 
equivalent hydrodynamic model, the Lagrangian-averaged Navier-Stokes (LANS) equations, also 
performed well in tests at -Ra ~ 300 (see |321] and references in D3I1 ). However, above Re « 3000 
(R\ ~ 800), it was shown that placing the filter width in the inertial range leads to contamination 
of the super- filter-scale properties (such as the spectra) for LANS. We refer here to this effect 
as the super- filter- scale bottleneck, which may be different in nature from the viscous bottleneck 
observed in some DNS of the Navier-Stokes equations. The contamination may be linked to the 
formation of spatial regions in the flow with no internal degrees of freedom (so-called "rigid bod- 
ies") J33I1 . which also correspond to the development of a secondary inertial range of the LANS 
equations at sub-filter scales. This super-filter-scale contamination provides an effective constraint 
on the filter size and, hence, on the available reduction of the total number of the (numerical) 
degrees of freedom (Ndof) needed to reproduce the large-scale dynamics of the flow at a given 
Reynolds number; a factor of ~ 10 can be achieved. This limitation is not apparent in low and 
moderate Reynolds number (resolution) simulations (e.g., 64 3 LANS compared with 256 3 DNS) 
as the scale separation is not enough for the above-mentioned phenomenon of contamination of 
small-scale spectra because of rigid body regions in the flow to appear. The bottleneck (and super- 
filter-scale contamination) was not studied as such but neither was it observed in 2D LAMHD for 



high Reynolds number 1127 



28 



29(1 . 3D LAMHD has only been tested at more moderate Reynolds 
number []30j] (see also Il34ll for a recent review). The aim of the present work is, thus, to determine 
if LAMHD in three space dimensions, for higher Reynolds number develops problems similar to 
that of LANS. Specifically, we test for the existence of spatial regions with no available internal 
degrees of freedom. We show in the following that LAMHD behaves better in this respect than 



LANS, and, thus, continues to appear as a promising model for MHD flows. 
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II. THE EQUATIONS OF MOTION 



We consider the incompressible MHD equations for a fluid with constant density, 

d t v + uxv=jxb - Vp + z/V 2 v 
9 ( b = Vx(vxb)+ r?V 2 b 

V-v = V-b = 0, (1) 

where v and b denote respectively the velocity and magnetic fields, p the pressure divided by 
the density, v the kinematic viscosity, and r/ the magnetic diffusivity. As is well known, in in- 
compressible MHD, Alfven waves will travel along a uniform background field, b . From linear 
perturbation analysis the dispersion relation between wavenumber, k, and frequency, uj, is 

(a; + i7]k 2 ) (uj + ivk 2 ) = k%l . (2) 

The wave speed is |b | and, assuming r] = u, the amplification factor is given by exp(— r]k 2 t). The 
ideal (77 = v — 0) quadratic invariants for MHD are in the L 2 norm. For example, the total energy 
is given by: 

Et = \ (NI2 + H&II2) = \^J d (\v\ 2 + |b| 2 ) d\. (3) 



The LAMHD equations Il25h are given by 



9 ( v + wxu = jxb- VII + zA7 2 v 
9 t b = Vx(uxb)+ r?V 2 b 
V-v = V- u = V- b = V- b = 0, (4) 

where u (b) denotes the filtered component of the velocity (magnetic) field and n the modified 
pressure. Filtering is accomplished by the application of a normalized convolution filter L : f \—> f 
where / is any scalar or vector field. By convention, we define u = v. LAMHD in the form given 
in Eqs. ([4]) is both computationally efficient and makes clear that Alfven's theorem is preserved by 
the model: the smoothed magnetic field is advected by the smoothed velocity. In the remainder of 
this paper, we take 77 = v (unit magnetic Prandtl number) and, thus, it is sufficient to introduce the 
same filtering for the velocity and magnetic fields in this case. This allows us to write LAMHD in 
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LES form, 



<9 t u + u>xu = jxb - VII + iA7 2 v - V ■ f 
9,b = Vx(uxb)+ r/V 2 b - V ■ f b 

V-v = V- u = V- b = V- b = 0. (5) 

We choose as our filter the inverse of a Helmholtz operator, L = Hr 1 = (1 — ct 2 V 2 ) -1 . 
Therefore, u = g a <g> v where g a is the Green's function for the Helmholtz operator, g a (r) = 
exp(— r / a) / (Ana 2 r) (i.e., the Yukawa potential), or in Fourier space, u(k) = v(k)/(l + a 2 k 2 ). 
The effective filter width is, thus, approximately a. With this choice, the Reynolds (turbulent) 
SGS stress tensor is given by 

f = a 2 (Vu ■ Vu T + Vu ■ Vu - Vu T • Vu - Vb ■ Vb T - Vb ■ Vb + Vb T ■ Vb) (6) 

and the divergence of the electromotive-force (emf) SGS stress tensor by 

V • f 6 = r/a 2 V 4 b. (7) 

In this form, the expression of the SGS tensors make explicit the fact that u = ±b Alfven waves 
are preserved even in the subgrid scales. These u = ±b waves travel along b (the smoothed and 
unsmoothed fields are identical for uniform b ) and the dispersion relation is 

(u + ivk 2 ) (u + iV k 2 (l + a 2 M k 2 )) = k 2 l t + a \ k l , (8) 

where a K and a M are the filter widths for the smoothing of the velocity and magnetic fields, 
respectively. For a = olk = % and r] = v (the case we study), the wave speed 
is given by b (l — {r\ka 2 k 2 /b ) 2 /8 + 0((r]ka 2 k 2 /b ) 6 )) , the strength of the smoothed back- 
ground magnetic field reduced by an order a 4 k 4 term. The amplification factor is given by 
exp (—T]k 2 t(l + a 2 k 2 /2)) for both u = —b waves traveling in the direction of b and u = b 
waves traveling anti-parallel to b . Finally, the ideal quadratic invariants for LAMHD are in the 
H^if) norm. For example, the total energy is given by a mixture of the smooth and rough fields, 
namely: 
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Ulrp 



|«) = II / ( u - a ! V 2 u) • u + (b - a 2 V 2 b) • b d> 
2 D J D 



x 



1 1 
2D 



v ■ u + b ■ b d x. 



(9) 



D 



We solve both sets of equations, Eqs. CD) a nd (J), for one specific instance of a decaying MHD 
flow, using a parallel pseudo spectral code [|35Ll36|l in a three-dimensional (3D) cube with periodic 
boundary conditions. The initial conditions for the velocity and magnetic fields are constructed 
from a superposition of three Beltrami (helical) ABC flows to which smaller-scale random fluc- 
tuations are added with initial kinetic and magnetic energy Ek = Em = 0.5, magnetic helicity 



Hm =< a • b 0.45 (b = V x a where a is the vector potential and the brackets denote vo 



ume 



38|] for 



average), and the initial co-alignment of the fields, (v • b) ( | v| |b| )~ 1 ps 10~ 4 (see Ill4 
details). A MHD-DNS with a resolution iV 3 = 1536 3 (i.e., 1536 grid points in real space in each 
direction) and rj — v — 2 • 10~ 4 is used as our high Reynolds number test case for the LAMHD 
model. The DNS computation is stopped when the growth of the total dissipation begins to enter 
the saturation phase (t = 3.7), at which time the Reynolds number based on the mechanical inte- 
gral scale is Re ~ 9200 and the Taylor Reynolds number w 1100. The MHD flow resulting from 
the initial conditions employed has previously been analyzed for its spectral properties and for the 
spatial structures it develops ||14ll37Ll38|l . In this paper, we perform a simulation with similar ini- 
tial conditions and parameters but now using LAMHD at a resolution of 512 3 grid points; we also 
perform for comparison purposes a Navier-Stokes LANS run with the same initial velocity field 
but with b = 0, on a grid of 512 3 points. In both cases, the filter width is a = 27r/18 (k a = 18) and 
is, thus, large enough to preclude any artifact of numerical resolution altering the results. Based on 



previous analyses |33j,|39J], we estimate k max /k" ps 2.4 (where k max is the maximum wavenumber 



resolved in the simulation, and k^ is the LAMHD dissipation scale) using computations conducted 
for rj = v = 6 • 10~ 4 with a Reynolds number of Re ~ 2200. However, the main point of using 
such a large filter is to test if LAMHD fails in the same way as LANS. We finally perform a LES 
simulation in a 256 3 grid using the LAMHD equations with the same viscosity and diffusivity as 
the 1536 3 DNS used for the comparison. In this way, we extend the Re ~ 9200 computation in 
time by a factor of 3. 
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III. RESULTS 



A. Spectral contamination in LANS for an ABC flow and its absence in the MHD case 

One of the main findings of our preceding work with LANS on the Navier-Stokes equations 
is that a k +1 scaling develops in the (kinetic) energy spectrum at sub-filter scales; this leads to 
a contamination of super-filter scales because of detailed energy conservation (per triadic inter- 
actions). This LANS k +1 spectrum (together with super- filter-scale spectral contamination) has 
only recently been recognized, in the case of one specific forcing function at large Reynolds num- 



ber [33], but such a spectral contamination has not yet been generally demonstrated (although 



theoretical arguments for the k +1 spectrum have been given in [33]). Thus, we first confirm its 
presence in a LANS simulation with the same viscosity and the (nearly) same initial conditions 
for the velocity field as for the MHD-DNS (and LAMHD runs) examined in this paper, and based 
on large-scale ABC flows with superimposed random noise at small scale. Due to the presence 
of random noise and considering the differences in resolution and the presence of a filter in the 
LAMHD runs, the initial conditions were not exactly reproduced, although the same procedure 
was used to generate them. In the present Navier-Stokes case, we find again what can be called 
an enhanced (super-filter-scale) bottleneck: the positive power-law spectral contamination of the 
kinetic energy spectrum E K {k) in the LANS run is observed for times after the peak of dissipa- 
tion (see dotted line, Fig. [Tk). The fitted spectrum is k +0 - 5 (note that k +1 requires the entire LANS 
spectrum to be resolved, and therefore has only been observed for much larger values of k max / k%). 

For the given parameters and initial conditions, we find the super- filter- scale bottleneck for 
LANS. However, when integrating the MHD equations with the Lagrangian model (dashed line, 
Fig. [It) with these same parameters, no such contamination is present. Note that the spectra for 
the DNS-MHD are shown at the time of peak dissipation, while the spectra for the Lagrangian- 
averaged models are for a slightly later time in order to allow for the possible formation of rigid 
bodies, which are known to be at the source of the spectral contamination close to the filter 
wavenumber in the Navier-Stokes case. For this reason, and due to the slight differences in initial 
conditions, we have chosen to plot spectra normalized to that of the DNS at k = 14 to emphasize 
the scaling. For most of the inertial range (also in an approximate sense below the filter width a) 
the scaling of E K (k) is reproduced by the LAMHD simulation. The sub-filter scaling for LAMHD 
is not as steep as MHD, but is not a positive scaling law. The agreement for E M {k) is remarkable. 
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FIG. 1: (a) Spectra of kinetic energy (normalized to DNS E K (U), see text) for 1536 3 MHD DNS (solid 
line), 512 3 LAMHD (dashed), and 512 3 LANS (dotted), in the latter case with b = at all times but 
otherwise identical conditions. For intermediate scales, k € [5,40], LAMHD reproduces the scaling of 
the DNS, the larger scales being affected by slight differences in initial conditions, see text. For k close 
to the filter scale (k € [k a /2, k a ]), a positive power law, A; 5 (gray line), is found for LANS, (b) Spectra 
of magnetic energy (normalized to DNS Em (14)) for the same runs: LAMHD reproduces the scaling of 
the DNS even beyond the filter wavenumber, k a = 18 as indicated by the vertical dashed line. LAMHD 
exhibits neither the positive power-law nor the super-filter-scale spectral contamination associated with high 
Reynolds number LANS modeling seen in (a). 



More importantly, neither positive-power-law spectra nor contamination of the super- filter- scale 
spectra are evidenced at all. 

B. The lack of rigid bodies in LAMHD in the large— a limit for unforced flows 

Evidence for the development of rigid bodies in LANS (which led to its limited use as a LES) 



has only been shown for I ^ a 1330. Since investigation of the large— a limit is not as computa- 
tionally demanding as the small—/ limit, it is interesting to look at this limit as a rough indication 
of what occurs for small a and smaller I. This approach has been employed both for the LANS 



Navier-Stokes case in two dimensions |40|] and in three dimensions p3fl . In such a case, the pur- 
pose is to examine the properties of the model itself, as opposed to trying to reproduce large-scale 
properties, the large-scale behavior being reduced to a very small span of wavenumbers. With this 
practice, the properties of the sub-filter-scales can be studied, to better understand the origin (or 
lack) of super-filter-scale contamination. We now use this limit to further explore the differences 
between LAMHD and LANS. We employ simulations for the two models with the same initial 
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conditions as before, with 77 = v — 5 • 10~ 5 (Re ~ 26, 000 at peak of dissipation for LAMHD), 
and a resolution of 256 3 grid points. Note that these dissipative coefficients are four times smaller 
than what was considered in the previous section since, for a fixed resolution, the achievable 
Reynolds number goes as a 2//3 . This follows for LANS from the predicted (and verified) degrees 



of freedom, Ng f oc a" 1 Re 3 / 2 [133 



39TI . The scaling of LAMHD may differ, but the same value of 



the viscosity is employed for the two models, regardless. 

For LANS, we observe the expected k +1 zero flux inertial range (see Fig. [2]) which is followed 
by a viscous (sub-filter-scale) bottleneck feature, k +1 - 5± - 2 , before the dissipative range proper. We 
conducted a second simulation with v = 10~ 4 and found a k 1A± ' 3 spectrum. This is analogous to 
results for DNS of the Navier-Stokes equations where only the viscous bottleneck is observed at 
moderate Reynolds number and is preceded by an inertial range only for higher Reynolds. These 
viscous bottlenecks may be different in nature from the (super-filter-scale) bottlenecks discussed 
before, which are not associated to the onset of the dissipative range but to the development of 
a secondary inertial range in LANS below the filtering length, and result in contamination of the 
large (resolved) scales when the LANS equations are used as an LES. Having confirmed that our 
analysis from the forced LANS case extends to the decaying LANS simulation, we now apply it to 
LAMHD. The large— a LAMHD spectra are given in Fig. |3] Notably, there is no positive-power- 
law spectrum. 



0.0100 x 




0.000 



1 
k 



FIG. 2: Spectrum of kinetic energy for a 256 3 grid with k a = 3 (v = 5- 10~ 5 ) LANS, b = (Navier-Stokes 
case). The fitted grey line, k +1A±A , agrees with the rigid-body hypothesis for the inertial range 13311 . This 
slope is followed by a steeper slope attributed to a bottleneck, with /c +15± - 2 . 



Predictions of energy spectra in the inertial range follow from the global scaling laws for third- 
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order structure functions for isotropic, homogeneous turbulence. Exact results for these structure 
functions have been found for incompressible MHD J41I and for LAMHD [29]. The latter are, in 
terms of both the smooth fields 7^ = u ± b and the rough fields = v ± b (where the z-fields 
are called the Elsasser variables): 



(10) 



where (.) denotes volume averaging, 5f(l) = /(x + 1) — /(x), and <5/||(l) = [f (x + 1) — f (x)] ■ 1. 
For sub-filter scales (I <C a), l 2 a 2 z ± and the scaling law becomes dimensionally zzz ~ el. 
This implies a sub-filter scale spectrum corresponding to the invariants E± = \\z ± \\ 2 t /2 for the 
ideal non-dissipative case. We then have E±(l)k ~ z ± z ± ~ (e±)' 



T a\2/3 a 2/3 or ^ equivalently, 



Em 



(ID 



as for LANS Q39Q. Recall that in the flux relation, Eq. (flOl) . e± stands for the energy transfer and 
dissipation rate of E±. Hence, the prediction, Eq. (fTTT) . for the spectra, E±(k), is, equivalently for 
E% = {\\uWl + ||6||q)/2 and for H% = ~ f D v • b (fx. The spectra shown in Fig. [3]for large— a 
LAMHD do not exclude, due to the large uncertainties of the fitted power laws, the predicted k~ l 
spectra. 




if 



10" 



io" 6 l , ! , , \, , , 

1 10 
fc 

FIG. 3: Spectra for a 256 3 grid with k a = 3 (7? = v = 5 ■ 10~ 5 ) LAMHD, Re « 26, 000: Total energy, 

Ex(k), (solid line) and cross helicity, Hc(k), (dashed). The fitted slopes, Er(k) ~ A;~°- 7± - 3 and Hc(k) ~ 
fc -0.5±.4 couW 

agree with either Kolmogorov or IK predictions for LAMHD (see text) at this level of 

uncertainty. 
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A spectral prediction for LAMHD can also be arrived at by dimensional analysis of the spec- 
trum which follows the scaling ideas originally due to Kraichnan Ii42l1 and which is developed for 
LANS in Ref. Il43n . Here, the energy dissipation rate, e± = dE±/dt, is related to the spectral 
energy density by 

' ' E%{k) (12) 



4 ~ 



where t fe is the turnover time for an eddy of size ~ k 1 . This turnover time is related to a "velocity," 



y±\2 



Z±Z±/(l + a 2 k 2 ) ~ kE«(k)/(l + a 2 k 2 ). Substitution 



Zf, (i.e., t fc ~ l/(kZ±)\ where (Z k 
into Eq. (PT2l) yields, 

E%(k) ~ (4) 2 / 3 Ar 5/3 (l + a 2 k 2 ) 1/3 



or, for «/c 3> 1, 



• A2/3 Jfe _1 . 



(13) 



(14) 



In the Iroshnikov-Kraichnan it 1 1 . 



1211 (hereafter, IK) phenomenology, Alfven waves (correspond- 



ing to either z T = 0) can only interact nonlinearly when they collide along field lines (along which 
they travel in opposite directions). The characteristic time for an Alfven wave is ~ (/cS ) _1 - If 
this is less than tfc, the effective transfer time i T is increased, i T ~ tl/U- Substitution of this new 
transfer time into Eq. (fT2l) yields, instead of Eq. (TT3T) 



El(k) ~ (45 ) 1/2 r 3 / 2 (i + a 2 k 2 Y 



/2 



(15) 



or, for «/c 3> 1, 



(16) 



The spectra shown in Fig. [3] for large— a LAMHD also agree with the IK predicted spectra, Eq. 
(fT6l) . In fact, the spectra more closely correspond to this prediction; this is consistent with the 
fact that, for this flow, an IK spectrum E(k) ~ k~ 3 ^ 2 is observed at large scale (followed by a 
weak turbulence anisotropic spectrum E(k±) ~ k~ 2 at small scale) [ 141]. Again, simulations at 
higher resolution are needed for a definite answer and the result may not be universal as shown 
for example in the context of reduced MHD dynamics due to the presence of a strong uniform 
magnetic field B 144|] or for MHD with a strong B 1131]. 

Another indication of the zero-flux regions in LANS is found by examining the spatial vari- 
ation of the cubed increments associated with the scaling laws 5uM)8ui(l)5vi(l) for LANS and 
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-10 -5 5 10 

FIG. 4: PDFs of cubed increments. The cubed increments when averaged are equal to flux times length, 
e a ■ I. Here I = 0.88a {a = 2vr/3). The dotted line is 5u\\(l)5ui(l)5vi(l) for LANS, solid for LAMHD 
5z^ 0-)Szf (l)Sz^ (1), and dashed for LAMHD 5zu~ {l)5z~ (l)Sz~ (1). More of the volume gives no contribu- 
tion to the flux for LANS than for LAMHD, indicating no rigid bodies in LAMHD. 

^(1)5^(1)^(1) for LAMHD (note that one can transform this relation into the u, v, b, b vari- 
ables). For a given length /, these cubed increments when averaged are related with the energy 
fluxes by Eq. (flOl) (the LANS relation and the hydrodynamic and MHD relations are contained in 
this expression in the corresponding limits). As a result of this correspondence, for brevity we will 
indicate cubed increments in the figures as the corresponding energy flux times the length used 
to compute the increments. This also allows us to identify regions with zero cubed increments as 
rigid bodies (a rigid rotation has zero longitudinal increments). Probability distribution functions 
(PDFs), see Fig. 0], indicate that LAMHD has a much smaller proportion of its volume, which 
could potentially be rigid bodies (i.e., frozen regions with no internal degrees of freedom (zero 
velocity increment), which therefore do not contribute to the energy flux). That is, more of the 
volume is contributing to the turbulent cascade. Snapshots for constructing the PDFs are taken 
from both a = 2n/3 Lagrangian-averaged models for times shortly after the peak of dissipation 
and when the LANS total dissipation is nearly equal to that of LAMHD. The strengths of the 
central peaks of the PDFs for large— a are another indication that LAMHD inherits none of the 
rigid-body or zero-fiux-region problems of LANS. 
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C. Why are spectral properties of LAMHD better than in the fluid case? 



Why does LAMHD not exhibit the same spectral contamination as LANS? One possible cause 
is the hyperdiffusivity term seen in the LES form for LAMHD, Eq. ([7]), whereas there is no 
hyperviscosity-like term in LANS. To test if this hyperdiffusion is responsible for the lack of 
spectral contamination in LAMHD, we removed the hyperdiffusion by setting f b = in Eqs. © 
or, equivalently, by substituting 7/V 2 b for 7/V 2 b in Eqs. ©. We then start the run from the same 
initial conditions but now with these new equations employing a = 27r/33 and v = rj — 2 • 10~ 4 
at a resolution of 384 3 (with hyperdiffusion, a smaller resolution of 256 3 is possible, see Section 
IIIIDI) . Note that such a modified LAMHD model is not expected to, nor found to, perform well as 
a SGS model; this numerical experiment is performed here only in order to assess the effect of the 
hyper-diffusive term introduced by the a modeling. We find that hyperdiffusion is not responsible 
for the lack of a k +1 spectral contamination in LAMHD (see Fig. [5]). 



0.1 000 

^ 0.0100 

0.001 
0.0001 

1 10 100 

k 

FIG. 5: Spectra for a 384 3 grid with k a = 33 obtained from the modified-LAMHD (see text) shortly after the 
maximum of dissipation: kinetic energy (solid) and magnetic energy (dashed); the LAMHD equations have 
been modified by removing the hyperdiffusive turbulent emf. Even without hyperdiffusivity, no positive 
power-law is found. Instead, fits (grey lines) for kinetic and magnetic energy spectra near the filtering 
length are k' 1 - 7 ^- 1 and fc~ 19±1 , respectively. 




Other possible causes for LAMHD not exhibiting the super- filter- scale bottleneck as does 
LANS are the actual physical differences between the two fluids that are modeled, Navier- 



Stokes and MHD. First, unlike incompressible Navier-Stokes, MHD supports oscillatory so 



(Alfven waves) which are linked to enhanced spectral nonlocality of energy transfer [a 



utions 



45D lead- 



ing to dynamic interactions between widely separated scales. For Navier-Stokes, the depletion of 
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energy transfer due to local interactions at some cutoff in wavenumber is believed to bring about 



the bottleneck effect 11461 . 



47 



48 



4911 . However, related to the spectrally nonlocal energy transfer 



via Alfven waves, MHD does not seem to exhibit a bottleneck in its spectra between the inertial 
and dissipative ranges 111411 . As LAMHD supports Alfven waves at all scales (and alters their dis- 
sipation and wave speed appreciably only for sub-filter scales), the same physics could be behind 
the lack of a super-filter-scale bottleneck in LAMHD. 

Another difference between the fluid and MHD cases is the geometry of the dissipative struc- 
tures: one finds vortex filaments for Navier-Stokes at high value of the vorticity, and current and 
vorticity sheets for MHD; sheets which are found to roll-up at high Reynolds number 13811 • It has 
been claimed that the development of helical filaments in the fluid case can lead to the depletion 



of nonlinearity and the quenching of local interactions 



5111 and, hence, to the viscous bot- 
tleneck. A similar energy transfer depletion may occur in LANS. In ^33\ evidence is presented 
that Taylor's frozen-in turbulence hypothesis applied to Lagrangian averages leads to the forma- 
tion of "rigid bodies" in the flow wherein there are no internal degrees of freedom and no transfer 
of energy to smaller scales (i.e. regions with e ~ 5u|/Z = as well as w x v = 0). These 
regions are likely related to the shorter, thicker vortex filaments formed and the suppression of 
vortex stretching dynamics as a is increased [52]. As MHD has spectrally non-local transfer (e.g., 
velocity at large scales does stretching of magnetic field lines at small scales) this leads to the 
break up of these rigid bodies in the LAMHD case and the breakup of the viscous bottleneck in 
the MHD case. The magnetic field interaction with the large scale velocity can re-enable transfer 
of energy to smaller scales of the velocity field. Indeed, defining the kinetic spectral transfer due 
to the Lorentz force as 



for LAMHD , and as 



T L (k) 



u fc • j x b dVt k 



v fc • j x b) dVt k 



(17) 



(18) 



for MHD, we see in Fig. [6] that the Lorentz force is removing large-scale kinetic energy and 
supplying small-scale kinetic energy; this effectively bypasses the formation of rigid bodies for 
LAMHD and the viscous bottleneck for MHD (note that Eqs. CLU) and ([18]) do not detail the 
scales at which magnetic energy is created or destroyed). 

This argument can also be recast in terms of Kelvin's circulation theorem. For Navier-Stokes, 
the circulation Y of the velocity v is conserved in the ideal case for barotropic flows. In ideal 
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FIG. 6: (Color online) Spectral transfer due to the Lorentz force, T L (for 1546 3 DNS) and T£ (for 512 3 
k a = 18 LAMHD) at a time just prior to the peak of dissipation. Positive Tl is shown as dash-dotted lines 
and negative Tl as dashed lines. Positive is shown as solid (green online) lines and negative T£ as 
dotted (green online) lines. LAMHD qualitatively reproduces the transfer of kinetic energy in MHD. 



MHD, this conservation is broken by the Lorentz force, 



-37 = -ji <J> v ■ dr = (j> j x h ■ dr, (19) 



dt dt 



where C is any material curve. As a result, while in ideal Navier-Stokes a material curve C defines 
the boundary of a vorticity tube with fixed strength, in MHD these structures are deformed and 
their vorticity content changed by the Lorentz force. A similar result follows for LAMHD and 
LANS, 

— = — (f> u ■ dr = ® j x b • <ir . (20) 
dt dt J c J c 

Breaking the conservation of circulation in this way can prevent the formation of a bottleneck. 
For example, for the fluid case in the Clark— a model (which differs from LANS only in the 
conservation of T), it was also found that no super-filter- scale bottleneck was present 115 311 . 



D. LES Application 

Having now shown that LAMHD does not suffer the same drawbacks with regards to energy 
spectra as LANS, we may turn our attention to a practical application. The purpose of a SGS model 
or LES is to make predictions about large Reynolds number flows at a reduced computational 
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FIG. 7: (Color online) Temporal evolution, r eddy ps 4.5, for 1536 3 DNS (solid, black), 256 3 k a = 33 
LAMHD (dashed, green online), 256 3 under-resolved "DNS" (dotted, red online), and 384 3 k a = 33 
nonhyperdiffusive-LAMHD (dash-dotted, blue online), (a) Time evolution of the energies: kinetic (lower 
curves), magnetic (middle curves) and total (upper curves), (b) Time evolution of total enstrophy, (j 2 + u 2 } 
((j 2 + u> • o>) for LAMHD and (j • j + us ■ u> ) for the nonhyperdiffusive case). Note that LAMHD gives 
a better agreement to the total dissipation rate up to the maximum time that the high resolution DNS is 
performed. Also note that the DNS equivalent to the LAMHD run presented here is not feasible on present- 
day computers at a reasonable cost. 



expense. From the scaling arguments in Refs. [|33l|39|], using simulations conducted at Re ps 2200, 
and assuming a k~ l scaling, we can estimate a = 1/33 for a 256 3 LAMHD-LES "prediction" of 
our 1536 3 MHD-DNS. Time evolution of the energies and the total enstrophy are shown in Fig. |7J 
for much later times than reasonably attainable with the MHD DNS with present-day computers. 
Also shown are results for solving the MHD equations, Eqs. (Q~|) with v — 2 ■ 10 -4 and a resolution 
of 256 3 : a so-called "unresolved DNS" and the non-hyperdiffusive modified-LAMHD from the 
previous section. Before the peak of dissipation, t ~ 4, the unresolved DNS gives a poorer 
prediction of the total dissipation and total energy which is then followed by a significantly larger 
and somewhat later peak of dissipation, at t ~ 5 than the resolved DNS and the LAMHD LES. 
The non-hyperdiffusive LAMHD is not expected to perform well as a SGS model and it is seen 
to be clearly under-dissipative. The ratio of magnetic to kinetic dissipation is ~ 1.5 for the DNS, 
ps 2.9 for LAMHD, pa 1.1 for the under-resolved DNS, and 1.4 for the non-hyperdiffusive model. 
Together with Fig. [7] (b) these ratios show that LAMHD achieves accurate total dissipation by 
an excess of magnetic dissipation and a reduction of kinetic dissipation (both at the small scales). 
This feature has already been depicted in Fig. 15 of Ref. |27|]. Compensated energy spectra 
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for the peak of dissipation (t £ [2.7, 3.7]) are shown in Fig. [8] For the under-resolved DNS, 
we observe the appearance of a tail at large wavenumbers with a k 2 spectrum as predicted using 
statistical mechanics arguments for truncated systems in the ideal (y = 0, rj = 0) case Il54ll . 
The under-resolved spectra are not significantly different from the resolved DNS, but note that a 
reliable and convincing determination of spectral indices, beyond visual inspection, does require 
high resolutions. Comparing now the resolved DNS and the LAMHD run, the quality of the 
spectra are similar for scales larger than a. Recall that differences at the largest scales, stem 
from the differences in initial conditions as stated in Section ITlI Al and from time evolution of the 
flow. Finally, noting that the computer saving here is 6 3 in memory and 6 4 in running time, we 
conclude that the LAMHD continues to behave satisfactorily, as already shown both in two space 



dimensions [27, 



28 



290 and in 3D OOh . in particular in the context of the dynamo problem of 



generation of magnetic fields by velocity gradients; thus, LAMHD may prove to be a useful tool 
in many astrophysical contexts where magnetic fields are dynamically important, such as in the 
solar and terrestrial environments, or in the interstellar and intergalactic media. 
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FIG. 8: Spectra compensated by k 3 / 2 for the kinetic (a) and magnetic (b) energies averaged over t € 
[2.7,3.7]; labels are as in Fig. |7] and the dashed vertical line indicates k a = 33. Note the k 2 tail at 
high wavenumber that is known to develop for under-resolved runs, a prediction stemming from statistical 
mechanics. 



We also computed a 512 3 LAMHD-LES (a = 1/ 85) which retains more of the small scales than 
the 256 3 LAMHD-LES while still yielding significant computational savings over the 1536 3 DNS. 
We compare this with the result for a = 1/18 (chosen not as a LES but to stress the model) in Fig. 
[9J The structure of sheets observed in MHD dissipative structures is preserved in the LAMHD 
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simulations, although current and vortex sheets become thicker in LAMHD as a result of the filter 
as a is increased. This is necessary to achieve reduced resolution computations. Note that these 
sheets are different in nature from the fat 'rigid bodies' observed in LANS, as the turbulent energy 
transfer to small scales is not quenched and there is no super-filter-scale bottleneck. 




FIG. 9: 2D cross sections of square current, j 2 , for 512 3 LAMHD-LES (a = 1/85) (Upper Left) and 
model-stress-case (a = 1/18) (Upper Right). MHD dissipative structures, sheets, are retained which 
become thicker as a is increased. (Lower Left) 256 3 LAMHD-LES (a = 1/33) and (Lower Right) 256 3 
unresolved DNS. For the unresolved run, current sheets are somewhat smeared out by numerical noise. 



IV. DISCUSSION 

In this paper, we have tested the LAMHD model against high Reynolds number direct nu- 
merical simulations (up to Reynolds numbers of ~ 9200) and in particular we have focused our 
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attention on the dynamics of small scales near the a cut-off. We find that the small-scale spec- 
trum presents no particular defect; specifically, we find that, unlike in the hydrodynamical case, 
the Lagrangian-averaged modeling for MHD exhibits, even at large Reynolds numbers, neither a 
positive-power-law spectrum nor any contamination of the super-filter-scale spectral properties. 
This difference between LANS and LAMHD is not due to the inclusion of a hyper-diffusive 
term in LAMHD that stems from the derivation of the model; rather, it stems from fundamental 
differences between hydrodynamics and MHD. Indeed, neither the (non-consistent) removal of 
hyperdiffusion from LAMHD nor the examination of scales much smaller than a gave any indi- 
cation of problems similar to those caused by the zero-flux regions found in computations using 
LANS. These regions limited the computational gains of using LANS as a LES in hydrodynamics 
to a factor of only 10 in computational degrees of freedom or 30 in computation time. LAMHD 
is not subject to the same limitations and, as we demonstrated, a gain of a factor of 200 in the 
number of degrees of freedom, or a factor of 1300 in computation time, obtains when comparing 
to the highest Reynolds number in turbulent MHD available today in a DNS. 

There are two obvious candidates to explain the lack of a (super-filter-scale) bottleneck effect 
in LAMHD: the enhanced (hyper-)diffusion in LAMHD compared with LANS, and physical dif- 
ferences between fluids and magneto-fluids, specifically, spectrally nonlocal transfer via Alfven 
waves and its associated breaking of the circulation conservation. The first candidate would elim- 
inate the super-filter-scale bottleneck by removing energy from the system and precluding the 
formation of a secondary range below the filtering scale a (note that this term becomes of the 
same order as the ordinary diffusion when I ~ a). Simulations of LAMHD performed without the 
hyper-diffusion term ruled out this scenario, as no super-filter bottleneck was found. 

The second candidate is the presence of the Lorentz force in MHD (and LAMHD) which breaks 
down the circulation conservation and provides the restoring force for Alfven waves. Both prop- 
erties were shown to be preserved by LAMHD. In Navier-Stokes, the development of helical 



filaments could quench local interactions [5 



3, 



51|1 depleting the energy transfer and leading to the 



viscous bottleneck. However, in MHD, the conservation of the circulation (dY/dt = in the 
absence of dissipation) is broken by the Lorentz force, which modifies Kelvin's theorem (see Eq. 
(fT9l)). The forcing term is associated with the Alfven waves, and represents the removal of circu- 
lation (and of kinetic energy) that is transfered to the magnetic field. Note that in Fourier space, 
the term scales as kE M (k) and is dominant compared to the dissipation in the inertial range. This 
term precludes the formation of rigid bodies, giving as a result a larger net flux towards smaller 



19 



scales and a resulting larger dissipation in MHD/LAMHD. This is illustrated in Fig. SI This sink 
of circulation may also be the cause of the lack of a viscous-scale bottleneck in MHD. In LANS it 



was shown [33, 



5311 that conservation of the circulation (except for viscosity) leads to the forma- 



tion of rigid bodies that fill a substantial volume of the fluid, and that in turn substantially decrease 
the energy flux to small scales, reduce dissipation, and create the super-filter scale bottleneck. In 
LAMHD, the destruction of sub-filter-scale rigid bodies by large scale magnetic field and shear 
results as the presence of a magnetic field permits the development of long-range interactions in 
spectral space fl&LZL|45|]- This can also explain why a— models for other non-local equations, or for 
problems that do not preserve the circulation provide good SGS models. As an example, the use of 
LANS in primitive equations ocean modeling gives satisfactory results, e.g. in its reproducing the 
Antarctic circumpolar current baroclinic instability that can be seen only at substantially higher 
resolutions when using direct numerical simulations 1^5 ]. 

Energy is dissipated in MHD flows through two different processes. Viscosity is responsible for 
the dissipation of mechanical energy, while Ohmic losses are responsible for dissipation of mag- 
netic energy. Mechanical and magnetic energy are not conserved separately, but rather coupled as 
illustrated by the existence of Alfven waves, which correspond to oscillations of the magnetofluid 
with the velocity field parallel or anti -parallel to the magnetic field, and associated to the inter- 
change of magnetic and kinetic energy. In MHD, it is believed that most of the total energy in the 
flow is finally dissipated (mediated by this interchange) through Ohmic losses, in a process that 
involves reconnection of magnetic field lines. This is supported by several simulations of MHD 



turbulence [5 



3 



570 and is consistent with phenomenology. While in hydrodynamics small scales 



are permeated by a myriad of vortex filaments, in MHD the dominant dissipative structures are cur- 
rent sheets, where strong gradients of the magnetic field and their associated strong currents lead 
to rapid Ohmic dissipation. Sub-grid models attempt to replace the physical processes of small- 
scale dissipation by processes that mimic the non-linear transfer of energy to smaller scales (where 
energy is in reality dissipated, but now in scales that are not resolved by the model). In traditional 
LES, this is done with enhanced turbulent viscosities. Note that the eddy viscosity is not obtained 
from the linear dissipative term (the term that describes the actual physical process responsible for 
the dissipation) but from the non-linear terms in the equations (the terms that describe the coupling 
between fields at different scales). The final goal is not to capture the dissipation processes, but to 
be able to preserve (with computational gains) the large scale dynamics. 

Lagrangian averaged models take a different (although related, see e.g., 12911 ) approach. Besides 
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adding (in some cases, as in the case of MHD) an enhanced viscosity, the non-linear terms are 
modified at small scales. This modification changes the time-scale of the energy cascade, and as 
a result changes the scaling law of the energy spectrum E(k) at sub-filter scales. This change 
leads to changes in the dissipation, as the dissipation is in the original equations proportional to 
k 2 E{k). The end result (an enhanced dissipation that is intended to mimic the transfer of energy 
to smaller scales in the unresolved scales) should be the same as in a traditional LES: gains in 
computing costs preserving as much information of the large scale flow as possible. As in the 
case of LES, the actual dissipation process is not as important as the fact that large-scale dynamics 
should be reproduced with minimal contamination by the sub-grid model. We believe the results 



presented here (and in earlier work 1127 . 



28, 



29. 



30, 



3111 ) show this is the case, and allow the 



use of the LAMHD equations as a subgrid model of MHD turbulence. However, considering 
the differences observed between LANS and LAMHD, we discuss the dissipation processes in 
LAMHD. Two mechanisms for dissipation can be identified in LAMHD: dissipation of mechanical 
energy through the viscosity, and dissipation of magnetic energy through (enhanced) Ohmic losses. 
From the equations, the total variation of energy goes as Q27D : dE/dt = —v (u> ■ u>) — r\ (j 2 ) and as 
a result the mechanical energy dissipation scales as k 2 Ev(k) while the magnetic energy dissipation 
scales as (1 + a 2 k 2 )k 2 E M (k). The extra k 2 factor in the latter gives more dissipation than in the 
LANS case. This excess of magnetic dissipation in LAMHD mimics, as previously mentioned, the 
dominant contribution to dissipation by Ohmic losses in MHD. This hyperdiffusion is required in 
the sub-filter scales to accurately model the total energy dissipated at the unresolved scales. This 
was demonstrated by our experiments with a modified LAMHD, where we (non-consistently) 
removed the hyperdiffusive term and found the resulting model to fail as a LES. 

Yet another way to understand the differences between LANS (for incompressible isotropic 
and homogeneous flows) and LAMHD is to consider the derivation of these models Il25n using 
the generalized Lagrangian-mean (GLM) formalism 115 811 . This form of Lagrangian averaging 
describes wave, mean-flow interactions. For the case of weak turbulence, where the nonlinear 
transfer is dominated by waves, GLM requires in principle no closure. As a result, GLM gives 
an exact closed theory for the evolution of the wave activity. On the other hand, when there are 
no waves (as in incompressible Navier-Stokes) or when eddies dominate the transfer, a closure is 
required. One possible closure assumes that fast fluctuations are just advected by the mean flow 
(basically, Taylor's frozen-in hypothesis for the small scale turbulent fluctuations) and leads to 
the several "a-models" that include LANS and LAMHD. In this context, it is not surprising for 
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subgrid models based on GLM to perform better in the presence of Alfven waves (for LAMHD) 
or Rossby and gravity waves (for the Lagrangian-averaged primitive equations [550). The more 
relevant the waves are to the dynamics, and to the non-linear coupling of modes in the system, the 
less relevant is the hypothesis behind the closure. Furthermore, the a-model equations can then 
be expected to be a better approximation to the problem at hand, that is, to be closer to an exact 
closure of the original system of equations. 

In the fluid case, the application of the "Taylor" closure that smaller- than- a scale fluctuations 
are swept along by the large-scale flow results in the fluctuations having greatly reduced inter- 
actions. This allows for a reduction in computational expense and leads to the super- filter- scale 
bottleneck by quenching spectrally non-local interactions. In the LAMHD case, the small-scale 
z + (z~) fluctuations are swept along by the large-scale z (z + ) flow. Small-scale fluctuations ad- 
vected by two different fields may now collide and nonlinearly interact. The second part of the 
model is the preferential hyperdiffusion of Alfven waves with wavelengths shorter than a. This 
damps rather than quenches nonlinear interactions among the small scales. This more gentle sup- 
pression of the transfer of energy to smaller scales reduces the numerical resolution requirements 
without forming a bottleneck. 

It was noted in [13011 when assessing the properties of LAMHD in the dynamo context that the 
overall temporal evolution was satisfactory, e.g. with a correct growth rate, although the growth 
of the magnetic seed field started slightly earlier in the LAMHD run than in the DNS. One can 
speculate as to whether this delay is linked to the super-bottleneck effect of LANS (which prevails 
when the magnetic field is negligible compared to the velocity, the two modeling approaches, 
LAMHD and LANS, being dynamically consistent). This point is left for future work; one could 
determine as well at what ratio of magnetic to kinetic energy the overshooting of spectra in LANS 
disappears for LAMHD. 

Also deserving of a separate study is to investigate the behavior of LAMHD when anisotropics 
that appear at small scales jl4 \ are present; this would be essential when a uniform magnetic 
field is imposed to the overall flow. The evaluation of the behavior of the model when computing 
spectra in the perpendicular and parallel directions (with respect to a quasi-uniform magnetic field, 
computed by locally averaging the field in a sphere of radius comparable to the integral scale) 
remains to be done but is somewhat time consuming. An analysis of the structures that develop in 
the highly turbulent LAMHD flow studied in the preceding section is also left for future work; of 
particular interest is the occurrence of Kelvin-Helmholtz like roll-up of current sheets as observed 
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at high resolution H14H ; however, the choice of the parameter a in the present paper was made 
on the basis of questioning the existence or lack thereof of a rigid-body high-wavenumber k +1 
spectrum and, thus, was not optimized for the study of the inertial range properties of the flow for 
which a much smaller value of the length a could be used. 

Finally, how far resolution can be reduced when using LAMHD as a LES for various statistics 
of interest will also require further detailed study. The present study shows that, to reproduce the 
super- filter- scale energy spectrum in three dimensions, gains by a factor of 1300 in computing 
time can be achieved. The need to reproduce higher order statistics can decrease these gains. 
As an example, in two-dimensional MHD, it was shown that gains when using LAMHD as a 
subgrid model depend for high order moments on the order that one wants to see to be accurately 
reproduced J29 I. 



Acknowledgments 

Computer time was provided by GWDG, NCAR, and the National Science Foundation Teras- 
cale Computing System at the Pittsburgh Supercomputing Center. The NSF Grant No. CMG- 
0327888 at NCAR supported this work in part and is gratefully acknowledged. PDM is a member 
of the Carrera del Investigador Cientifico of CONICET. The anonymous referees are gratefully 
acknowledged for improving the clarity of the discussion of our results. 



[1] C. Meneveau and J. Katz. Scale-Invariance and Turbulence Models for Large-Eddy Simulation. An- 
nual Review of Fluid Mechanics, 32:1-32, 2000. 

[2] A. Pouquet, U. Frisch, and J. Leorat. Strong MHD helical turbulence and the nonlinear dynamo effect. 
Journal of Fluid Mechanics, 77:321-354, 1976. 

[3] A. Yoshizawa. Subgrid modeling for magnetohydrodynamic turbulent shear flows. Physics of Fluids, 
30:1089-1095, 1987. 

[4] J.-P. Chollet and M. Lesieur. Parameterization of Small Scales of Three-Dimensional Isotropic Tur- 
bulence Utilizing Spectral Closures. Journal of Atmospheric Sciences, 38:2747-2757, 1981. 

[5] J. Baerenzung, H. Politano, Y. Ponty, and A. Pouquet. Spectral modeling of magnetohydrodynamic 
turbulent flows. Phys. Rev. E, 78(2):026310-+, August 2008. 



23 



[6] A. Alexakis, P. D. Mininni, and A. Pouquet. Shell-to-shell energy transfer in magnetohydrodynamics. 
I. Steady state turbulence. Phys. Rev. E, 72(4):046301-n, 2005. 

[7] P. Mininni, A. Alexakis, and A. Pouquet. Shell-to-shell energy transfer in magnetohydrodynamics. II. 
Kinematic dynamo. Phys. Rev. E, 72(4):046302-i-, 2005. 

[8] Y. Zhou, O. Schilling, and S. Ghosh. Subgrid scale and backscatter model for magnetohydrodynamic 
turbulence based on closure theory: Theoretical formulation. Phys. Rev. E, 66(2):026309— n, 2002. 

[9] S. Galtier, S. V. Nazarenko, A. C. Newell, and A. Pouquet. A weak turbulence theory for incompress- 
ible magnetohydrodynamics. Journal of Plasma Physics, 63:447-488, 2000. 
[10] P. Goldreich and S. Sridhar. Toward a theory of interstellar turbulence. 2: Strong alfvenic turbulence. 
ApJ, 438:763-775, 1995. 

[11] P. S. Iroshnikov. Turbulence of a Conducting Fluid in a Strong Magnetic Field. Soviet Astronomy, 
7:566^, 1964. 

[12] R. H. Kraichnan. Inertial-range spectrum of hydromagnetic turbulence. Physics of Fluids, 8:1385- 
1387, 1965. 

[13] J. Mason, F. Cattaneo, and S. Boldyrev. Numerical measurements of the spectrum in magnetohydro- 
dynamic turbulence. Phys. Rev. E, 77(3):036403-+, 2008. 

[14] P. D. Mininni and A. Pouquet. Energy Spectra Stemming from Interactions of Alfven Waves and 
Turbulent Eddies. Physical Review Letters, 99(25):254502-i-, 2007. 

[15] Annick Pouquet, Pablo Mininni, David Montgomery, and Alexandras Alexakis. Dynamics of the 
Small Scales in Magnetohydrodynamic Turbulence . In Yukio Kaneda, editor, Proceedings of the IU- 
TAM Symposium on Computational Physics and New Perspectives in Turbulence, Nagoya University, 
Nagoya, Japan, September, 11-14, 2006, volume 4, pages 305-312. Springer- Verlag, 2008. 

[16] O. Agullo, W.-C. Miiller, B. Knaepen, and D. Carati. Large eddy simulation of decaying magnetohy- 
drodynamic turbulence with dynamic subgrid-modeling. Physics of Plasmas, 8:3502-3505, 2001. 

[17] M. L. Theobald, P. A. Fox, and S. Sofia. A subgrid-scale resistivity for magnetohydrodynamics. 
Physics of Plasmas, 1:3016-3032, 1994. 

[18] N. E. L. Haugen and A. Brandenburg. Hydrodynamic and hydromagnetic energy spectra from large 
eddy simulations. Physics of Fluids, 18:5106—1-, 2006. 

[19] W.-C. Miiller and D. Carati. Dynamic gradient-diffusion subgrid models for incompressible magne- 
tohydrodynamic turbulence. Physics of Plasmas, 9:824-834, 2002. 

[20] D. W. Longcope and R. N. Sudan. Renormalization group analysis of reduced magnetohydrodynamics 



with application to subgrid modeling. Physics of Fluids B, 3:1945-1962, 1991. 

[21] B. Knaepen and P. Moin. Large-eddy simulation of conductive flows at low magnetic Reynolds num- 
ber. Physics of Fluids, 16:1255-1-, 2004. 

[22] Y. Ponty, P. D. Mininni, D. C. Montgomery, J.-F. Pinton, H. Politano, and A. Pouquet. Nu- 
merical Study of Dynamo Action at Low Magnetic Prandtl Numbers. Physical Review Letters, 
94(16): 164502-+, 2005. 

[23] Y. Ponty, H. Politano, and J.-F. Pinton. Simulation of Induction at Low Magnetic Prandtl Number. 

Physical Review Letters, 92(14): 144503-+, 2004. 
[24] D. D. Holm. Averaged Lagrangians and the mean effects of fluctuations in ideal fluid dynamics. 

Physica D Nonlinear Phenomena, 170:253-286, 2002. 
[25] D. D. Holm. Lagrangian averages, averaged Lagrangians, and the mean effects of fluctuations in fluid 

dynamics. Chaos, 12:518-530, 2002. 
[26] D. C. Montgomery and A. Pouquet. An alternative interpretation for the Holm "alpha model". Physics 

of Fluids, 14(9):3365-3366, 2002. 
[27] P. D. Mininni, D. C. Montgomery, and A. G. Pouquet. A numerical study of the alpha model for 

two-dimensional magnetohydrodynamic turbulent flows. Physics of Fluids, 17:5112—1-, 2005. 
[28] J. Pietarila Graham, P. D. Mininni, and A. Pouquet. Cancellation exponent and multifractal structure 

in two-dimensional magnetohydrodynamics: Direct numerical simulations and Lagrangian averaged 

modeling. Phys. Rev. E, 72(4):045301(R)-+, 2005. 
[29] J. Pietarila Graham, D. D. Holm, P. Mininni, and A. Pouquet. Inertial range scaling, Karman- 

Howarth theorem, and intermittency for forced and decaying Lagrangian averaged magnetohydro- 
dynamic equations in two dimensions. Physics of Fluids, 18:045106, 2006. 
[30] P. D. Mininni, D. C. Montgomery, and A. Pouquet. Numerical solutions of the three-dimensional 

magnetohydrodynamic a model. Phys. Rev. E, 71(4):046304 — 1-, 2005. 
[31] P. D. Mininni. Turbulent magnetic dynamo excitation at low magnetic Prandtl number. Physics of 

Plasmas, 13:056502-+, 2006. 
[32] Alexey Cheskidov, Darryl D. Holm, Eric Olson, and Edriss S. Titi. On a Leray— a model of turbulence. 

Proceedings of the Royal Society of London, A46L629-649, 2005. 
[33] J. Pietarila Graham, D. Holm, P. Mininni, and A. Pouquet. Highly turbulent solutions of 

the Lagrangian-averaged Navier-Stokes alpha model and their large-eddy-simulation potential. 

Phys. Rev. E, 76:0563 10-+, 2007. 



[34] P.D. Mininni, A. Pouquet, and P. Sullivan. Two examples from geophysical and astrophysical turbu- 
lence on modeling disparate scale interactions. In Roger Temam and Joe Tribbia, editors, Summer 
school on mathematics in geophysics. Springer- Verlag, 2008. to appear. 

[35] D. O. Gomez, P. D. Mininni, and P. Dmitruk. MHD simulations and astrophysical applications. Ad- 
vances in Space Research, 35:899-907, 2005. 

[36] D. O. Gomez, P. D. Mininni, and P. Dmitruk. Parallel Simulations in Turbulent MHD. Physica Scripta 
Volume T, 116:123-127, 2005. 

[37] W. H. Matthaeus, A. Pouquet, P. D. Mininni, P. Dmitruk, and B. Breech. Rapid Alignment of Velocity 
and Magnetic Field in Magnetohydrodynamic Turbulence. Physical Review Letters, 100(8):085003— v, 
2008. 

[38] P. D. Mininni, A. G. Pouquet, and D. C. Montgomery. Small-Scale Structures in Three-Dimensional 
Magnetohydrodynamic Turbulence. Physical Review Letters, 97(24): 244503— v, 2006. 

[39] C. Foias, D. D. Holm, and E. S. Titi. The Navier-Stokes-alpha model of fluid turbulence. Physica D 
Nonlinear Phenomena, 152-153:505-519, May 2001. 

[40] E. Lunasin, S. Kurien, M. Taylor, and E. Titi. A study of the Navier-Stokes-alpha model for two- 
dimensional turbulence. ArXiv Physics e-prints, 2007. 

[41] Helene Politano and Annick Pouquet. Dynamical length scales for turbulent magnetized flows. Geo- 
physical Research Letters, 25(3):273-276, 1998. 

[42] R. H. Kraichnan. Inertial Ranges in Two-Dimensional Turbulence. Physics of Fluids, 10:1417-1423, 
1967. 

[43] C. Cao, D. D. Holm, and E. S. Titi. On the Clark a model of turbulence: global regularity and 

long-time dynamics. Journal of Turbulence, 6:19-+, 2005. 
[44] P. Dmitruk, D. O. Gomez, and W. H. Matthaeus. Energy spectrum of turbulent fluctuations in boundary 

driven reduced magnetohydrodynamics. Physics of Plasmas, 10:3584-3591, 2003. 
[45] A. Alexakis. Nonlocal Phenomenology for Anisotropic Magnetohydrodynamic Turbulence. ApJ, 

667:L93-L96, 2007. 

[46] J. R. Herring, D. Schertzer, M. Lesieur, G. R. Newman, J. P. Chollet, and M. Larcheveque. A compara- 
tive assessment of spectral closures as applied to passive scalar diffusion. Journal of Fluid Mechanics, 
124:411^37, 1982. 

[47] D. Lohse and A. Miiller-Groeling. Bottleneck Effects in Turbulence: Scaling Phenomena in r versus 
p Space. Physical Review Letters, 74:1747-1750, 1995. 



[48] D. O. Martinez, S. Chen, G. D. Doolen, R. H. Kraichnan, L.-P. Wang, and Y. Zhou. Energy spectrum 
in the dissipation range of fluid turbulence. Journal of Plasma Physics, 57:195-201, 1997. 

[49] P. D. Mininni, A. Alexakis, and A. Pouquet. Large-scale flow effects, energy transfer, and self- 
similarity on turbulence. Phys. Rev. E, 74(1):016303-+, 2006. 

[50] H. K. Moffatt and A. Tsinober. Helicity in laminar and turbulent flow. Annual Review of Fluid 
Mechanics, 24:281-312, 1992. 

[51] Arkady Tsinober. An Informal Introduction to Turbulence. Kluwer Academic Publishers, Dordrecht, 
2001. 

[52] S. Chen, D. D. Holm, L. G. Margolin, and R. Zhang. Direct numerical simulations of the Navier- 

Stokes alpha model. Physica D Nonlinear Phenomena, 133:66-83, 1999. 
[53] J. Pietarila Graham, D. D. Holm, P. D. Mininni, and A. Pouquet. Three regularization models of the 

Navier-Stokes equations. Physics of Fluids, 20(3):035107-n, 2008. 
[54] U. Frisch, A. Pouquet, J. Leorat, and A. Mazure. Possibility of an inverse cascade of magnetic helicity 

in magnetohydrodynamic turbulence. Journal of Fluid Mechanics, 68:769-778, 1975. 
[55] M. W. Hecht, D. D. Holm, M. R. Petersen, and B. A. Wingate. Implementation of the LANS-a 

turbulence model in a primitive equation ocean model. Journal of Computational Physics, 227:5691- 

5716, 2008. 

[56] N. E. L. Haugen, A. Brandenburg, and W. Dobler. Is Nonhelical Hydromagnetic Turbulence Peaked 

at Small Scales? ApJ, 597:L141-L144, 2003. 
[57] P. D. Mininni. Inverse cascades and a effect at a low magnetic Prandtl number. Phys. Rev. E, 

76(2):026316-f, 2007. 

[58] D. G. Andrews and M. E. Mclntyre. An exact theory of nonlinear waves on a Lagrangian-mean flow. 
J. Fluid Mech., 89:609-646, 1978. 



27 



